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Traditionally, trends of universal turbulence statistics are presented versus Rx, which is the 
Reynolds number based on Taylor's scale A and the root-mean-squared (rms) velocity Urms- A 
and Urms, and hence R\, do not have the attribute of universality. The ratio of rms fluid-particle 
acceleration to rms viscous acceleration, _Ra, is an alternative to R\ that has the advantage of being 
determined by the small scales of turbulence. This ratio has the following attributes: _Ra is a 
Reynolds number, it is composed of statistics of the small scales of turbulence, can be evaluated 
with single- wire hot-wire anemometry, and like R\, can be partially evaluated by means of flow 
similarity. For isotropic turbulence the relationship between i?a and R\ is given. Graphs of 
velocity derivative flatness measured in anisotropic turbulence driven by counter rotating blades 
have different appearances depending on whether R^ or Rx is the abscissa. 



I. INTRODUCTION 

Reynolds sought, from the Navier-Stokes equation, "the dependence of the character of motion on a relation 
between the dimensional properties and the external circumstances of motion." Assuming that the motion depends 
on a single velocity scale U and length scale c, Reynolds found that the accelerations are of two distinct types 
and thereby deduced that the relevant solution of the Navier-Stokes equation "would show the birth of eddies to 
depend on some definite value of cpU j [i^'' j|] where p is the mass density of the fluid and [l is the coefficient of 
viscosity. Reynolds performed exhaustive experiments that demonstrated his deduction, as well as experiments on 
the stabilization of fluctuating flow. [Q He discovered the sudden onset of flow instability The Navier-Stokes 
equation is a = 9u/9t-t-u- Vu — — Vp-|-i^V^u, where -p is pressure divided by p, i/ = pj p is kinematic viscosity, u is the 
velocity vector, and a is the acceleration. Batchelor [|] discussed (in his Sec. 4.7) the interpretation of the Reynolds 
number as a measure of a relative to the viscous term vV^w.. He noted that the balance of the Navier-Stokes equation 
can also be parameterized in terms of the relative magnitudes of Vp and v'^^w. The latter parameterization does 
not technically lead to a Reynolds number, but it will be shown that the two parameterizations become equivalent at 
large Reynolds numbers. We therefore introduce the two ratios: 

i?a^(a.a)i/V(^Mv'u)-(V2u))'/' and R^, ^ {^p ■Vv)"'' I {v^ {V^ry) ■ (1) 

where (•) denotes an average. -Ra is a Reynolds number in Batchelor's concept, and i?vp is a parameterization of the 
relative magnitudes of Vp and z^V^u. 

To paraphrase Nelkin's ||^ description of Reynolds number scaling: if two turbulent flows have the same 
geometry and the same Reynolds number, then their statistics, when appropriately scaled, should be equal. A 
statistic of the small scales of turbulence is an average of quantities that contain only products of differences, such as 
two-point velocity difference or derivatives of velocity. Universality of the small scales of turbulence is the hypothesis 
that statistics of the small scales, when appropriately scaled, should become equal as Reynolds number increases |^ 
[Q; that is, the flow geometry becomes negligible in the limit that the Reynolds number is infinite. Discovering the 
appropriate scaling that results in universality is the topic of a vast amount of research |^] Q and will not be pursued 
here. The relevance of universality to real turbulent flows is discussed by Nelkin jsj and Sreenivasan and Antonia [Q. 

The Reynolds number based on the root-mean-square (rms) of the longitudinal-velocity component Wrms = 

1/2 I ^'^ 

(u\) and Taylor's length scale A is R\ = Urms^jv, where v is kinematic viscosity, and A = Urms/ \ (dui/dxi) ) 

Here, ui and xi are the components of velocity and spatial coordinate in the direction of the 1-axis. For decades, Rx 
has been used as the abscissa for presenting statistics that are believed to be universal aspects of small-scale turbulence 

(such as velocity derivative statistics normalized by powers of (^{dui / dxif'^) . The observed trends as R\ increases 
are an often-sought quantification of scaling universality. R\ has the advantage of being easily measured because 
it requires only measurement of ui (which yields dui/dxi by means of Taylor's hypothesis); that measurement can 
be obtained with a single hot-wire anemometer. Alternatively, flow similarity can be used to estimate the energy 
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dissipation rate e, and by substituting the local-isotropy relationship that e = 15iy (^(dui/dxi)'^J, R\ can be obtained 

from R\ = u^ms/ (^'^/IS)^''^. Because R\ depends on Urms, it depends on large-scale geometry of the flow. Nelkin 
discussed the nonuniversal attributes of R\. As a result of the nonuniversality of Rx, statistics of the small scales, 
e.g., normalized derivative moments, when graphed with Rx on the abscissa, can have different curves corresponding 
to dissimilar flows. 



II. ALTERNATIVE 



In addition to graphing such statistics with Rx on the abscissa, it would seem advantageous to use a quantity 
on the abscissa that is solely a property of the small scales of turbulence. That advantage has long been recognized. 
lH l@l Here, we consider the alternative i?a defined in (]^), and determine how it can be measured with an 
instrument no more complex than a single- wire hot-wire anemometer. Because the intended application is to statistical 
characteristics of the small scales, it is appropriate to simplify i?a and Rx on the basis of local isotropy. Indeed, local 
isotropy is a precondition for universahty. Q On this basis, (a • a) — (Vp • Vp) + iv'^ (V^u) • (V^u)) |^ f^; in 
which case 

i?a = + R%^ and i?vp = ^JRI - 1. (2) 
Because (Vp • Vp) » (v^ (^^'^) ' (^^")) ^^^S^ Reynolds numbers [||, (||) becomes i?a — R\i-p- 



A. Menu for the variance of viscous acceleration 



The correlation of viscous acceleration is ||] 

(r) ^ ly^lu^^l.v!:) = 



-V^V^A, (r) 



(3) 



where prime denotes evaluation at a point x'; x and x' are independent variables; r = x — x'; is the Laplacian 
operator in r-space; the right-most expression in (|^) is obtained on the basis of local homogeneity. Let e denote 
the energy dissipation rate, and (r) and I?ijfe (r) denote the second- and third-order velocity structure functions. 
The Navier-Stokes equation and local isotropy give [0 dtDiAv) -f dr^Dijk (r) -I- (4/3) e = 2vV^Dij (r); applying the 



(4) 



operator — to that equation and comparison with gives, 

v., (r) = - {v/A) Vl [dtD,, (r) + 9., (r)] 



Derivative operators are abbreviated; e.g., dt — d/dt, = d'^/dr^, etc. Summation is implied by repeated indexes. 
Performing the contraction of (||) and such that the terms become functions of r and — > 9^ + (2/'') dr, we have, 
on the basis of local isotropy. 



(r) 



dt Df3p (r) 



24, 



rd^ + lld:^ + —df. Dii{r) 



-dr 



Du (r) + 



ri^lll (r) + -d- 



(5) 



Power-series expansion of the structure functions followed by differentiation and taking the limit r — > gives various 
formulas for the viscous acceleration variance: 

y,. (0) = (V^u) . (V^u)) = U.'{{dlu,y) = 35.^ {{dlu,] 



(6) 

where = d^/dx^. The term — (1/2) dte, which vanishes for local stationarity, is included above, but was ne- 
glected by Ref. ||]. The enstrophy generation rate can be written as: (—35/2) (^{dxiUi)^'^ = {ujiUjjSij) — 

(—4/3) {sijSjkSki)', comparing this with the right-most expression in (^ gives expressions that can be evaluated by 
DNS or multi-wire anemometers Ol 11121: 
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V^^ (0) = (V^u) • (V^u)) = -^fte + u {ujiujjs.j) = -\dte - {s^jSjkSki) , (7) 
where uji is vorticity vector and Sij is rate-of-strain tensor. The velocity-derivative skewness is 

S EE (^{d,,u,f) I {{d,,u,f)"^ = {{d,,u,f) I l^ElX^vf^ , (8) 

such that gives 

(V^u) • (V^u)) - -\^dte + 0.3£3/2j.-i/2 \s\ . (9) 

The variation of l^l with R\ is known. Herring and Kerr ||l3| show \S\ increasing from 0.074 at R\ — 0.46 to 
become constant at 0.5 for R\ > 20 (see Ref. [Q), whereas at R\ > 400 Antonia et al. show j^l increasing from 0.5 

as \S\ ~ 0.8i?5I'"'^^. If e and Urms are known from measurements or from flow similarity, then R\ = u^^g/ {ev/lbf'^'^ 
and hence 15*1 is known such that (^) can determine (y'^ fV^u) ■ (V^u)). The term — (1/2)9(6 can be neglected 
if it is not evaluated. The DNS data of Herring and Kerr [jisl |0 suggest that \S\ ~ R\/h for R\ < 1, then 
(V^u) • (V^u)) ~ 0.06e3/2z.-i/2i?^ for R^ < 1. 



B. Menu for the variance of pressure-gradient acceleration 

Poison's equation, local homogeneity and local isotropy, but no other approximations, result in fl^ ] 

oo 

(Vp • Vp) = 4 J [i^ini (r) + (r) - 6Dii,^p (r)] dr, (10) 



where (r), Daaaaii"), and I?ii/3/3 (r) are components of the fourth-order velocity structure-function tensor, 

which is defined by Tiij^i (r) = (^{ui — u'^) (uj — u^) {uk — it'j.) (u; — mJ)); the 1-axis is parallel to the separation vector 
r; a and /? denote the Cartesian axes perpendicular to the 1-axis. Thus, a and /? are 2 or 3; equally valid options 
under local isotropy are a — P or a ^ l3. 

There is enough cancellation between the positive 

and negative parts of the integrand, i.e., between [-Dim {"r) + Daaaa i^)] and —r~^6Diif3/3 (r), to make evalu- 

oc 

ation of the integral / [-Dim (r) + Daaaa (r) ~ QDupp (r)] dr difficult by means of experimental or DNS data 



£1 il 111- Hill and Wilczak jl| argued that the ratio =/ r^^ [D^^j^^ (r) + Daaaa {r) - GDnpp (r)] dr/ J 



oo 

r^^Diiii (r) dr = (Vp • Vp) /4 J r^^Dim (r) dr is a universal constant at high Reynolds numbers. Universality 



oc 

of is equivalent to the assertion that (Vp ■ Wp) scales with / r^'^-Dim (r) dr at high Reynolds numbers. Hill 



and Wilczak jl^ pointed out that the utility of determining is that the pressure-gradient variance can then be 
measured with a single-wire hot-wire anemometer by means of 

oc 

{Vp-Vp)^AH^ j r-^D^i^i{r)dr. (11) 



Asymptotic expressions for (Vp • Vp) derived from ([ll| ) are [|l8| 

(Vp • Vp) ~ i.lH^e^/'^v-^l'^F^ '^'^ ~ 3.9iJ^e3/2j,-i/2^o.25 f^^. > 4qq^ (-j^2) 

where F = {ipxiUi)'^^ j (^(S^jUi)^^ , and the right-most expression is obtained from the data for F of Antonia et al. 
[0, and, for low Reynolds numbers [|l^ 

(Vp • Vp) ~ 0.11e3/2j,-i/2^^ foj. < 20. (13) 
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Using DNS data, the preferable evaluation of is via = (Vp • Vp) / 4 / r '^-Diiii ir) dr so as to avoid 

- J-T3i 

the statistical uncertainty caused by the cancellations within the integrand of (|lC). Vedula and Yeung evaluated 
using DNS data with R\ < 230 and obtained the variation from ~ 0.55 at R\ = 20 to an approach to a 
constant value of ~ 0.65 in the range 80 < R\ < 230. For iJy = 0.65, ( |l^ ) agrees quantitatively with the DNS 
data in Table 1 of Gotoh and Fukayama for Rx > 387, and ( p^ ) is a good approximation for Rx > 200. Since 
the asymptotic formula ( ]T^ ) is thereby verified with = 0.65, the range of validity of = 0.65 seems to extend to 
the highest Rx for which Antonia et al. obtained F, namely Rx ~ 10^, but it might extend to Rx — oo. On that 
basis, we have the increase from the low-Reynolds-number asymptote — 0.36 as Rx — > 0, to ~ 0.55 at 

Rx = 20, to - 0-65 in the range Rx > 80. 

Thus, the menu for calculating (Vp • Vp) is (p^), or (11) with the empirically known values oi as a function 
of Rx, or ([l^)-(|l3|) with the DNS data of Vedula and Yeung [|l9| and Gotoh and Fukayama |2^ used in the range 
20 > i?A > 400. 



III. ISOTROPIC TURBULENCE 



For isotropic turbulence there is a one-to-one relationship between i?a and Ry, it is shown in Fig. 1 with i?vp 
included. Fig. 1 was obtained from (||) by use of (|l^)-([l3|) and (||) and the discussions following those equations. 

That is: for Rx < 1, Rvp ^ {l2Rx/ib \S\f''^ ~ (12/7)^/^ ~ 1.3; for 1 < i?A < 20, i?vp ^ {l2Rx/ib \S\f''^ with data 

of Kerr Q for \S\] for 20 < i?A < 400, i?vp - 2.6 ((Vp • Vp) / e^/^v'^''^f''^ with data tabulated by Vedula and Yeung 

and Gotoh and Fukayama [|o|; for Rx > 400, i?vp - [(3.1i/,^e3/2j,-i/2^o.79) / i^q_^^3/2^^i/2 \s\)]^^^ ~ 3.3i?^o^; 

where the above mentioned Rx variation of F and \S\ of Ref. |7| was used; finally, R^ = (l + R^pf^ . In Fig. 1, 
i?Vp — 3.37?° °'' is extended to Rx = 100 to show, by comparison with the data of Vedula and Yeung and Gotoh and 
Fukayama, that it is a good approximation for Rx > 200. From i?A = 1 to lO'' in Fig. 1, i?vp changes by one decade 
and i?a less so. Fig. 1 is based on nearly isotropic data only for Rx < 400; for Rx > 400 the data used for F in 
( p^ ) and S in (||) is that of Antonia et al. [Q which consists, in part, of atmospheric surface layer data at Rx > 2000. 
The assumption is that for _Ra ^ 400 the turbulence is locally isotropic and that the relationships of F and S to Rx 
measured by Antonia et al. |0] do not differ significantly from what they would be for isotropic turbulence. Some 
support for the assumption is that the asymptote i?vp — 3.3i?" '^^ in Fig. 1 agrees with the nearly isotropic data to 
Rx — 200. However, futher confirmation must await DNS and experiments on nearly isotropic turbulence at higher 
Reynolds numbers than have been attained to date. 



A. Relationship to recent data 

The DNS data of Vedula and Yeung [|I9| , Gotoh and Rogallo E2| , and Gotoh and Fukayama ||2^] are in the 
range of Rx where \S\ ~ 0.5. Since their DNS are steady state, (|) gives {i^^ (V^u) • (V^u)) ~ 0.15£^/^i^~^/^. 
In their Fig. 1, Vedula and Yeung ||l^ show a ratio that they call C, which equals R^p, as well as a quantity 
= (a ■ a) / (3e3/2j,-i/2) 0.05RI, where (i^^ (V^u) • (V^u)) ~ 0.15e3/2jy-i/2 ^as used; they graphed both C and 

a^^ versus Rx- Similarly, Gotoh and Rogallo |2|] and Gotoh and Fukayama Q show F\/p — 3ao^^ ~ 0.15i?a versus 
Rx in their Figs. 1 and 2. Therefore, the above-mentioned graphs show R^p and i?^ as Rx varies. Reversing the 
role of ordinates and abscissas in their graphs and in Fig. 1, the graphs show Rx for nearly isotropic turbulence as 
the universal Reynolds number i?a varies. 



IV. ANISOTROPIC TURBULENCE 

For anisotropic turbulence, the value of (uf ) depends, by definition of anisotropy, on the direction of r. Thus, 
Rx does also. The appearance of a statistic graphed versus Rx can change depending on which velocity component 
is used to calculate Rx- For many flow geometries, both i?vp and i?a increase as Rx increases. It is appropriate to 
change perspective: The nonuniversal abscissa Rx typically, but not necessarily, increases as the universal abscissas 
R\/p and _Ra increase. 

The turbulent flow in a cylinder driven by counter rotating blades has several modes of large-scale structure 
p3| ||2^ and is anisotropic at large scales, but local isotropy is approached (25) at high Reynolds numbers. Belin 
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et al. 1^ measured F in the flow between counter-rotating blades; their F graphed with R\ on the abscissa has a 
maximum and minimum, shown here in Fig. 2a. That behaviour of F has been called controversial |^7[] because it 
has not been previously observed, and was not observed by Pearson and Antonia who obtained F in the range 
of R\ where Belin et al. observe the maximum and minimum. The Belin et al. data for F is thus a good test of 
whether or not i?a produces a different organization of a small-scale statistic than does R\. Only the single lowest 
Reynolds-number datum of Belin et al. is below R\ = 200. Figure 1 shows that (|l^ ) is accurate for R\ = 200. From 

the menus we therefore select and ( p^ and use = 0.65; then i?vp — (2.0i^°'^®/0.3 jS*!)^^^. Thus their data for 
F and S are used to calculate i?vp, then i?a is obtained from (^, although i?a — Rvp in this case. In Fig. 2a the 
data of Belin et al. for F versus R\ is shown using distinct symbols in several ranges of R\. Figure 2b shows their 
F versus i?a using the same symbols for the same data points. The maximum and minimum in Fig. 2a is absent in 

1 /2 

Fig. 2b, and the symbols are significantly rearranged in Fig. 2b relative to Fig. 2a. Figure 2b is not just (f"'^^) 
vs. F because that would be a straight line. The data is scattered in Fig. 2b because of the scatter in S. Their 
5- values at nearly the same value of R\ vary by as much as 28%. The important point is shown by comparing Fig. 
2a with 2b: The use of R^ on the abscissa rather than R\ can give a significantly different distribution of the data 
of a small-scale statistic. 



V. SUMMARY AND COMMENT 

As defined in (^, i?a is a Reynolds number; it is the ratio of rms fluid-particle acceleration to rms viscous 
acceleration; it is composed of statistics of the small scales of turbulence; it can be used as a universal abscissa for 
judging the universality of turbulence statistics. Although i?vp is not strictly a Reynolds number, it can also be used 
as a universal abscissa. At high Reynolds numbers, i?vp — ^a- Rvp can be evaluated with single- wire hot-wire 
anemometry, and by flow similarity; use of (||) then determines i?a- If one chooses to measure S for the purpose of 
obtaining the viscous acceleration variance from then accurate values of S require many samples and high spatial 
resolution. This does not require an enormous data set because the data can be stored as a histogram having one 
bin for each of the digitizer's distinct outputs with a bin for digitizer overflow and one for underflow. 

Models of the small-scale statistics of turbulence should be expressed in terms of universal attributes instead 
of in terms of R\. For example, in Table II of Belin et al. |2^], the model by PuUin and Saffman |2^ is in good 
agreement with data when judged in terms of power laws between derivative moments, but it is in relatively poor 
agreement with data when judged in terms of power laws between normalized derivative moments and R\. R\ can 
be specific to the flow geometry. 
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Fig. 1 Relation of i?a (upper symbols) and i?vp (lower symbols) to R\ for isotropic turbulence. Data of Refs. 
iQ, [|l9|, produce the pluses, asterisks, and triangles, respectively. Solid lines are the low- and high- Reynolds 
number asymptotes. 

Fig. 2 a) The velocity-derivative flatness of Belin et al. |^ versus R\ with different symbols for several ranges of 
Rx; b) the same data and symbols for each point, but versus i?a- 
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